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The scenario of baryogenesis through leptogenesis is reviewed in models involving heavy 
Majorana neutrinos. The various mechanisms of CP violation occurring in the out- 
of-equilibrium lepton-number-violating decays of heavy Majorana neutrinos are studied 
within a rcsummation approach to unstable- particle mixing. It is explicitly demonstrated 
how the resummation approach preserves crucial field-theoretic properties such as uni- 
tarity and CPT invariance. Predictions of representative scenarios are presented after 
solving numerically the Boltzmann equations describing the thermodynamic evolution 
of the Universe. The phcnomenological consequences of loop effects of heavy Majorana 
neutrinos on low-energy observables, such as lepton-flavour and/or lepton-number non- 
conservation in r and Z-boson decays and electron electric dipolc moment, are discussed. 

1. Introduction 

Present astronomical observations related to abundances of the light elements 4 He 
and 4 He, the content of protons versus antiprotons in cosmic rays, etc, lead to 
the conclusion 1,2 that, before the nucleosynthesis epoch, the Universe must have 
possessed an excess in the baryon number B which is expressed by the small baryon- 
to-photon ratio of number densities 

^ = (4-7) x HT 10 . (1.1) 

This baryonic asymmetry, uab = tib — ti§ ~ ?i_b, should have survived until today 
if there had been no processes that violate the B number and/or modify the number 
density of photons n 7 . Sakharov, assuming that the Universe was created initially 
in a -B-conserving symmetric state, was able to derive three necessary conditions 3 
to explain the baryon asymmetry in the Universe (BAU): 

(i) There must be -B-violating interactions in nature, so that a net B number can 
in principle be generated. 
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(ii) The i3-violating interactions should violate the discrete symmetries of charge 
conjugation (C) and that resulting from the combined action of charge and 
parity transformations (CP). In this way, an excess in baryons over antibaryons, 
AB, is produced. 

(iii) The B- and CP-violating processes must be out of thermal equilibrium, 4,5 
namely they should have an interaction rate smaller than the expansion rate 
of the Universe. This last requirement ensures that the produced AB is not 
washed out by the inverse processes. 

Grand unified theories (GUT's) can in principle contain all the above neces- 
sary ingredients for baryogenesis. 6 In such theories, out-of-equilibrium B- and CP- 
violating decays of super-heavy bosons with masses near to the grand unification 
scale Mx ~ 10 15 GeV can produce the BAU. However, this solution to the BAU has 
its own problems. The main difficulty is the generic feature that minimal GUT's 
predict very small CP violation, since it occurs at very high orders in perturba- 
tion theory. This problem may be avoided by augmenting GUT's with extra Higgs 
representations. 1 Also, GUT's must comply with limits obtained by experiments on 
the stability of the proton. Such experiments put tight constraints on the masses of 
the GUT bosons mediating B violation and their couplings to the matter. Another 
severe limitation to scenarios of baryogenesis arises from the anomalous B + L- 
violating processes, also known as sphalerons, 7 ' 8 ' 9 which are in thermal equilibrium 
for temperatures 10 ' 11 200 <~ T <~ 10 12 GeV. Unlike B + L, sphalerons preserve the 
quantum number B—L. Therefore, any primordial BAU generated at the GUT scale 
should not rely on B + L- violating operators, which imposes a further non-trivial 
constraint on unified theories. In that vein, Kuzmin, Rubakov and Shaposhnikov 9 
suggested that the same anomalous B + L- violating electroweak interactions may 
produce the observed excess in B during a first-order electroweak phase transition. 
Such a mechanism crucially depends on the Higgs-boson mass M ff , 12 ' 13 and the 
experimental fact Mh > 80 GeV practically rules out this scenario of electroweak 
baryogenesis. 14 Therefore, baryogenesis provides the strongest indication against 
the completeness of the SM, as well as poses limits on its possible new-physics 
extensions. 

Among the many baryogenesis scenarios invoked in the literature, the most 
attractive one is due to Fukugita and Yanagida 15 , and our emphasis will be put 
on their scenario in this review article. In such a scenario, out-of-equilibrium L- 
violating decays of heavy Majorana neutrinos Ni, with masses mjv 4 3> T c , produce 
an excess in the lepton number L which is converted into the desired excess in B by 
means of B + L-violating sphaleron interactions, which arc in thermal equilibrium 
above the critical temperature T c . Over the last years, many authors have discussed 
such a scenario, also known as baryogenesis through leptogenesis. 16 ' 17 ' 18 ' 19 ' 20 ' 21 
However, we should remark that isosinglet heavy neutrinos are not theoretically 
compelling in order to create an excess in the L number. Recently, Ma and Sarkar 22 
suggested a leptogenesis scenario based on a generalized Higgs-triplet model, 23 ' 24 
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Fig. 1. One-loop (a) self-energy and (b) vertex graphs in heavy Majorana neutrino decays. 



where the leptonic asymmetry is generated by out-of-cquilibrium CP-violating de- 
cays of heavy doubly charged Higgs triplets into charged leptons. However, such 
alternatives seem to face the known gravitino problem if they are to be embedded 
in a supersymmetric theory. 25 The charged Higgs triplets or their supersymmetric 
partners may interact strongly with gravitinos and produce them in large abun- 
dances. The slow decay rate of gravitinos during the nucleosynthesis epoch distorts 
the abundances of the light elements at a level inconsistent with present observa- 
tions. 

Mechanisms that enhance CP violation play a decisive role in baryogenesis. 
Using the terminology known from the K°K° system, 26 one may distinguish the 
following two cases: 

(i) CP violation originating from the interference between the tree-level decay 
amplitude and the absorptive part of the one-loop vertex. Such a mechanism 
is usually called e'-type CP violation. 15 ' 16 ' 17 ' 18 

(ii) CP violation induced by the interference of the tree-level graph and the ab- 
sorptive part of a one-loop self-energy transition. This mechanism is termed 
e-type CP violation. 27 ' 28 ' 9 ' 19 ' 20 ' 21 

As can be seen from Fig. 1, both of the above two mechanisms of CP violation 
are present 19,20,21 in the usual leptogenesis scenario 15 of heavy Majorana neutrino 
decays. CP violation of the e' type was extensively discussed in the literature. 
15,16,17,18 jf a vj Yukawa couplings of the Higgs fields to N and the ordinary lep- 
ton isodoublets are of comparable order, 16 ' 18 then baryogenesis through the e'-type 
mechanism requires very heavy Majorana neutrinos with masses of order 10 7 -10 8 
GeV. Such a high mass bound may be lifted if a strong hierarchy for Yukawa cou- 
plings and to^ is assumed. 16,17 However, without the latter assumption, 18 one 
obtains e' < 10 -15 for 1 TeV, and hence very heavy neutrinos are needed to 

account for the BAU. 

Recently, e-type CP violation and its implications for the BAU has received 
much attention. 20,21 In particular, it has been observed 20 ' 21 that CP violation can 
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be considerably enhanced through the mixing of two nearly degenerate heavy Ma- 
jorana neutrinos. Such an analysis cannot be performed in the conventional field- 
theoretic framework, since finite-order perturbation theory breaks down in the limit 
of degenerate particles. To be specific, the wave-function amplitude that describes 
the CP-asymmetric mixing of two heavy Majorana neutrinos, Ni and N2, say, is 
inversely proportional to the mass splitting mjvi — mjv 2 , and it becomes singular if 
degeneracy is exact. Solutions to this problem have been based on the Weisskopf 
and Wigner (WW) 29 approximation, 20 and the resummation approach. 30 ' 21 Both 
approaches lead to similar conclusions concerning the resonant enhancement of CP 
violation. Here, we shall follow the latter method, as the discussion of many cru- 
cial field-theoretic issues, such as renormalization, CPT invariance and unitarity, is 
conceptually more intuitive in this framework. 

To describe the dynamics of CP violation through mixing of two unstable par- 
ticles, one is compelled to rely on resummation approaches, which treat unstable 
particles in a consistent way. 31 ' 32 In fact, to any finite order in perturbation the- 
ory, physical amplitudes reflect the local gauge symmetry, respect unitarity, are 
invariant under the renormalization group, and satisfy the equivalence theorem. 
All of the above properties should also be present after resummation. Unfortu- 
nately, resummation methods often end up violating one or more of them. The 
reason is that subtle cancellations are distorted when certain parts of the ampli- 
tude arc rcsummcd to all orders in perturbation theory, whereas others, carrying 
important physical information, are only considered to a finite order. In this con- 
text, a novel diagrammatic resummation approach has been developed, 33 which is 
based on the pinch technique (PT) 34 and devoid of the above pathologies. In the 
PT resummation approach, basic field-theoretic requirements, such as analyticity, 
unitarity, gauge invariance and renormalizability, 33 are naturally satisfied. Apart 
from the great phcnomcnological importance of such a resummation formalism for 
the proper definition of the mass and the width of unstable particles, such as the W, 
the Z boson and the Higgs boson, 33 ' 35 this formalism may also be extended to the 
case of mixing between two intermediate resonant states in scattering processes 36 ' 30 
retaining all the required field-theoretic properties mentioned above. 

The afore-mentioned resummation formalism has been proved to be very suc- 
cessful in describing resonant transitions taking place in collider experiments. These 
are situations where the unstable particles are produced by given asymptotic states, 
e + e~, say, and their subsequent decay is observed by detecting some other asymp- 
totic states in the final state, e.g. e + e~ — > Z* — > However, in an expanding 
Universe, the unstable particles may undergo a huge number of collisions before 
they eventually decay. Each of these collisions contributes a Coulomb phase shift, 
and hence the mixed heavy particles are practically uncorrelated when they decay. 
To some extent, this thermodynamic phenomenon may be described by Boltzmann 
equations. 4 In this context, a related formalism for decays has been developed, 21 
which effectively takes into account decoherence phenomena in the mixing and sub- 
sequent decay of heavy particles, namely heavy Majorana neutrinos in our case. 
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Specifically, it is shown that £-type CP violation can even be of order unity. 21 This 
is in agreement with earlier studies on resonant CP violation through mixing in 
scatterings involving top quarks, supersymmetric quarks or Higgs particles in the 
intermediate state. 31 ' 36 > 30 Finally, wc must remark that alternative formulations 
of Boltzmann equations already exist in the recent literature 37 but they are ex- 
pected not to alter drastically the existing conclusions 20 ' 21 as far as the resonant 
phenomenon of CP violation is concerned. 

The organization of the review article is as follows: in Section 2 we briefly 
review the basic theoretical background concerning the B + L anomaly in the Stan- 
dard Model (SM), and the effect of sphaleron processes on the chemical potentials 
of SM particles. These considerations lead to a relation between the generated 
leptonic asymmetry and the observed baryonic asymmetry induced by sphaleron 
interactions. In Section 3 we discuss theories that naturally include heavy Majo- 
rana neutrinos. For illustration, we consider a minimal model with two isosinglet 
neutrinos and demonstrate how CP and L can simultaneously be violated in this 
model. In Section 4 we address the issue of renormalizability of the minimal iso- 
singlet neutrino model. Section 5 discusses in detail the resummation approach and 
its effective extension to describe incoherent decays of heavy unstable fermions. In 
Section 6 we apply the effective approach to the decays of heavy Majorana neutri- 
nos. In Section 7 we explicitly demonstrate how the resummation approach satisfies 
unitarity. In Section 8, we solve numerically the Boltzmann equations for repre- 
sentative leptogenesis scenarios, and give numerical estimates and comparisons for 
the BAU generated through e- and/or e'-type CP violation. Furthermore, we esti- 
mate the impact of finite-temperature effects on the resonant phenomenon of CP 
violation. Heavy Majorana neutrinos may also have important phenomenological 
implications for low-energy observables, as they can give rise to a non-vanishing 
electric dipole moment (EDM) of the electron at two loops or induce L-violating 
decays of the Z boson and the r lepton. These new-physics effects are detailed in 
Section 9. Section 10 summarizes our conclusions. 

2. B + L anomaly and sphaleron processes 

In the SM, the B and L numbers are only conserved in the classical action. After 
quantization, however, both baryonic and leptonic currents are violated by triangle 
anomalies, i.e., 

M = « = ^ (-a w W^ a W^+a Y Y^ v ), (2.1) 

where Np is the number of flavours, and a w = 3^/(4^), a Y = g 2 /(47r), are the 
SXJ(2) L and U(l)y fine-structure constants, respectively. Similarly, W pv , Y^ v 
are their respective field-strength tensors, and the antisymmetric tensors W^ v — 
\znv\ P W Xp , = \e llv \ p Y xp are their associated duals. Furthermore, baryonic 
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and leptonic currents are denned as 

J b = \ £^W (2.2) 

q,a> 

Jl = + (2-3) 

where q, I and denote quarks, charged leptons and neutrinos, respectively, and 
the index a indicates the colour degrees of freedom of the quarks. Since Eq. (2.1) 
also holds for individual lepton families, the actual anomaly-free charges are 

\b - L e , l -B - L M) Ib - L T . (2.4) 

It is then obvious that B + L symmetry is anomalously broken at the quantum level. 

The different gauge field configurations are characterized by different Chcrn- 
Simons numbers ncs- The CS numbers label the infinitely many degenerate vacua 
of the system. The variation of B + L number due to a quantum tunnelling from 
one vacuum state into another is given by 

A(B + L) = 2N F ^ / d A x W^' a W" = 2N F An CS • (2.5) 

o7T J 

At zero temperature, 't Hooft 7 estimated the probability of B-violating processes, 

and found them to be extremely suppressed by a factor exp(— 47rncs /ce w ) w exp(— 150ncs) 

relative to the ^-conserving ones with ncs = 0. 

The situation changes drastically at finite temperatures. The effect of non- 
trivial topological instanton-type solutions, termed sphalerons, 8 is amplified at high 
temperatures, thereby enhancing also the rate of the B- violating processes. To be 
precise, sphaleron interactions are in thermal equilibrium for temperatures in the 
interval 

100 GcV ~ T ~ 10 12 GeV. (2.6) 
Sphalerons may be thought of as the creation out of the vacuum of a state 

n i= i,jv F (uLd L d L v L )i . (2.7) 

Since these interactions violate B + L, any primordial baryonic asymmetry B should 
have a significant component in B — L or in the charges stated in Eq. (2.4), which 
are preserved by sphalerons, whereas any B + L component will be washed out. 
Decays of heavy Majorana neutrinos produce an excess in L, which can naively be 
written as a sum of an excess in \{B + L) and in ^(B — L). Sphalerons will then 
erase the B + L component but preserve the B — L, so one expects that about half 
of the leptonic asymmetry L will be converted into the baryonic asymmetry B, and 
also be preserved as B — L asymmetry. As we will see below, a more careful analysis 
based on chemical potentials leads to the conclusion that sphalerons approximately 
convert one-third of an initial leptonic asymmetry L into the observed baryonic 
asymmetry B. 
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For illustration, wc shall assume that all SM particles are almost massless at 
temperatures above the critical temperature T c . Actually, they have thermal masses 
but, to leading order, these are small and may be neglected. The number density 
of a particle a is given by 

f d 3 p a 1 

n " 9a J 2E a {2irf exp[(25 Q - » a )/T\ ± 1 ' 

where g a counts the internal degrees of freedom of a, p a and E a = (\p a \ 2 + m^ 1 ! 2 
are the three-momenta and the energy of the particle, respectively. The plus sign 
in Eq. (2.8) is for particles obeying the Fcrmi-Dirac statistics and the minus for 
particles governed by the Bose-Einstein statistics. The chemical potential for anti- 
particles, e.g. that of a, is opposite to that of the particles, i.e. p a — —p & - The later 
relation is valid if particles and antiparticlcs have interaction rates with photons or 
other gauge particles much higher than the expansion rate of the Universe. This is 
almost the case for all SM particles. However, this is not generally true for non-SM 
particles, such as isosinglet or right-handed neutrinos, which do not have any tree- 
level coupling to the W- and Z- bosons; their couplings are suppressed by loops and 
small Yukawa couplings. Under these assumptions, the number-density asymmetry 
of a SM particle a versus its antiparticle a is easily estimated by 

"At. = n a - n & w T 3 (J^) . (2.9) 

We shall now turn to an analysis of chemical potentials in the SM. Since FCNC 
interactions are sufficiently fast, we assign the same chemical potential for all differ- 
ent families of up and down quarks, i.e. (p UL , pd L , Muhi In contrast to quarks, 
individual leptons possess different chemical potentials, i.e. {pi L , p UlL , pi R ), where 
I = e, /i, r. Furthermore, the chemical potential of all neutral gauge bosons, such 
as gluons, photons, and Z bosons, vanish, and pw is the chemical potential of the 
W~ boson. Finally, the components of the Higgs doublet [x~, <p° = (H — ix°)/V2] 
have chemical potentials (po,p-). Many chemical potentials can be eliminated by 
means of chemical equilibrium reactions in the SM. More explicitly, we have 

W~ <-» 
W~ <-» 
w- <-» 

0° <-> 
0° «-» 
0° «-» 

As independent parameters, we consider p u — n UL , p — J2i ^il = J2i Ml > Mo 
and p,w- In the SM with Np families and Nh Higgs doublets, the baryon and 
lepton number B and L as well as the electric charge Q and hypercharge Qz may 
be expressed in terms of these quantities, as follows: 



u L + d L , 
vil + II , 

x- + 0° , 

U_ L + U R , 

d,L + d Rl 
II + Ir, 



Pw 
Pw 
Pw 
Po 
Pa 
Po 



— PvtL + Wl) 

P- + Po, 

-Ml + Vdni 



(2.10) 



B = 4N F p u + 2N F p w , 
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L = 3/U + 2Np^w — NpHo > 

Q = 2^^-2^ + 2(2^ + ^)^0-2(2^ + 2 + ^)^, 
Q 3 = -{2N F + A + 2N H )^ W . (2.11) 

Furthermore, the sphaleron interactions in Eq. (2.7) give rise to the additional 
relation 

N F (3nu + 2n w ) + /i = 0. (2.12) 

Above the electroweak phase transition, both charges Q and Q3 are conserved, i.e. 
(Q) = (Q3) = 0. Thus, we have: [iw = 0, [i = —3N F /j, u , and fi = —8N F [i u / (4N F + 
2Nh). Using these relationships among the chemical potentials, it is not difficult 
to obtain 38 ' 39 

8N F + 4Nh 

BIT > T c ) = F — (B-L) . (2.13) 

1 ' 22N F + 13N H V ' V ' 

From Eq. (2.13), one concludes that almost independently of the number of genera- 
tions and Higgs doublets, roughly one-third of the initial B — L and/or L asymmetry 
will be reprocessed by sphalerons into an asymmetry in B. This amount of B asym- 
metry persists even after the electroweak phase transition. 

3. Models with heavy Majorana neutrinos 

GUT's such as SO(10) 40,41 or E 6 42 models naturally predict heavy Majorana 
neutrinos. These theories also contain several other particles, e.g. leptoquarks, 
additional charged and neutral gauge bosons, etc., which may have significant in- 
teractions with heavy Majorana neutrinos and so affect the number density of heavy 
neutrinos. To avoid excessive complication, we shall assume that these new parti- 
cles are much heavier than the lightest heavy Majorana neutrino, and are therefore 
expected to decouple sufficiently fast from the process of leptogenesis. 

As already mentioned, SO(10) 40,41 and/or E 6 42 models may naturally accom- 
modate heavy Majorana neutrinos. Specifically, SO(10) models can break down to 
the SM gauge group in the following schematic way: 

SO(10) -» G 422 =SU(4)p S ®SU(2) i? ®SU(2) L 

-> G 3221 = SU(3) C ® SV(2) R ® SU(2) L ® U(1) (B _ L) 

-> SM = G 32 i =SU(3)®SU(2) L ®U(l)y, (3.1) 

where the subscript PS characterizes the Pati-Salam gauge group. 43 The spinor 
representation of SO (10) is 16-dimensional and its decomposition under G422 reads 

G422 : 16 -» (4,1,2) (4,2,1). (3.2) 

Evidently, SO(10) contains the left-right-symmetric gauge group SU(2)i{®SU(2)£,® 
U(1)(b_l), which necessitates the presence of right-handed neutral leptons. In 
this scenario, there can exist several Higgs-boson representations that may cause a 
breaking of the groups G422 and G3221 down to the SM gauge group G321. 41,44 
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Eg theories 42 may also have a breaking pattern related to SO(10) theories. In 
fact, the 27 spinor representation decomposes into 16 © 10 1 under SO(10). This 
leads to four singlet neutrinos per SM family: one neutrino as isodoublct member 
in 16, two neutrinos as isodoublet members in 10, and one singlet neutrino in 1. 
In these models, it is argued that depending on the representation of the Eq Higgs 
multiplets, two of the four isosinglets can have Majorana masses of few TeV, 42 
whereas the other two may be very heavy with masses of the order of the unification 
scale. 

We shall now discuss a generic subgroup that may be derived from SO(10) and/or 
Eg models. This generic subgroup may be realized in the usual SM augmented by 
a number hr of right-handed neutrinos um, with i = 1,2, ... ,ur. As we have 
discussed above, in Eg theories the active isosinglet neutrinos may be more than 
three. In SO(10) models, left-right symmetry is more naturally realized with one 
right-handed neutrino per family (for interesting alternatives, see Ref. 41 ). For the 
sake of generality, we shall keep the number of iso-singlet neutrinos arbitrary. For 
definitcness, the quark sector of the minimal model has the SM form, while the 
leptonic sector consists of the fields: 

Ir , v Rl , 

where I = e, fi, t. At T > T c ~ v, the vacuum expectation value (VEV) of the SM 
Higgs doublet <I> at temperature T (with v = v(0j) vanishes, v(T) — 0. At these high 
temperatures, all SM particles including Higgs fields are massless; they only acquire 
thermal masses. However, one may have Majorana masses in the Lagrangian given 

by 

~£m = 2 E ( p Ri M iJ^ + vmM»*v c Rj ). (3.3) 

where M v is an ur x n^-dimensional symmetric matrix, which is in general complex. 
In Eq. (3.3), the superscript C denotes the operation of charge conjugation, which 
acts on the four-component chiral spinors ipL and tpR as follows: 

{^l) C = PrC^ t , (Va) C = P L C4, T , (3.4) 

where Pl{r) = [1 ~ (+)7s]/2 is the chirality projection operator. The mass matrix 
M v can be diagonalized by means of a unitary transformation 

U T M V U = M v , (3.5) 

where U is an tir x n^-dimensional unitary matrix and the diagonal matrix M u 
contains the Ur heavy Majorana masses. Then, the respective tir mass eigenstates 
Ni are related to the flavour states vr% through 

ur n R 

VRl = Pr^U^Nj, v c m = Pl^U^Nj. (3.6) 

J=l 3=1 
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In the mass basis of heavy Majorana neutrinos, the Yukawa sector governing the 
interactions of the heavy neutrinos with the Higgs and lepton isodoublets is given 

£y = -EE ^ ^ ( {H ~S )/V * ) N * + H - c - ( 3 - 7 ) 

1=1 3 = 1 \ X / 

At high temperatures, the CP-even Higgs field H, the CP-odd Higgs scalar x° 
and the charged Higgs scalars given in (3.7) are massless. In the low-T limit 
T <C T c , the field H becomes the massive SM Higgs boson, whereas x° an d X ± arc 
the would-be Goldstone bosons eaten by the longitudinal degrees of freedom of the 
gauge bosons Z and W^, respectively. In the calculations in Section 6, we shall 
include the Mh dependence in the CP asymmetries. 

Let us now consider a simple one-generation model where the standard fermionic 
content is extended by adding two isosinglet neutrinos, e.g. vr and Sl- Then, the 
most general Yukawa sector that preserves lepton number reads 

-C = \{S L , (v R f ) ( ^ ^ ) ( {S ^ C ) + h R {p L , I l )*vr + H.c, (3.8) 

where $ = i<72<& is the isospin conjugate Higgs doublet and a-i is the usual Pauli 
matrix. The kinematic parameters M and fiR may in general be complex but their 
phases are not physical. One can make both real by appropriate redefinitions of the 
fermionic fields, i.e. 

L T L = ( VL ,l L )^e^L T L , VR ^e^v R , S L -> e» B S L . (3.9) 

One choice could be: <f>R — 0, <f>s — ar g(-^)7 an d <f>L = arg(Zi^). Retaining the 
L-conserving structure of the isosinglet mass matrix, such scenarios require a non- 
trivial mixing among the generations to describe CP violation. 45 Furthermore, such 
scenarios do not produce the necessary leptonic asymmetry for baryogencsis; how- 
ever, see Ref. 46 for an interesting variant based on individual lepton-flavour viola- 
tion. 

In order to break both L and CP symmetries of the Lagrangian in Eq. (3.8), 
one must consider at least two extensions in the model: 

(i) The inclusion of two complex L-violating Majorana masses /-IrDrVr and 

(ii) The addition of the L-violating coupling Hl {pL, Il)&(Sl) C and one //-violating 
mass parameter, e.g. hrvrv^. 

The two models are related by a unitary rotation and are therefore equivalent. The 
necessary conditions for CP invariance in these two scenarios may be found to be 



(i) |M 2 Im(M*Wfl) = 0, 

(ii) lm(h L h* R n R M*) = 0. 



(3.10) 
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It is now interesting to remark that fiL and Hr can be much smaller than M 
within E 6 scenarios. 42 These parameters may be induced after integrating out high- 
dimensional operators involving ultra-heavy non-active neutrinos. One may think 
that the lepton number is somehow violated, at the GUT or Planck scale Mx, by 
these additional non-active isosinglct fields, and it is communicated to the active 
isosinglet sector where M <C Mx- In this way, one can naturally obtain a see-saw- 
like relation for the sizes of and [1r, i.e. 

M 2 M 2 
HL, M* ~ jj^ or — , (3.11) 

where Ms ~ 10~ 3 Mx could be some intermediate see-saw scale. In such generic 
mass models, the heavy Majorana neutrinos Ni and 7V 2 have a very small mass 
splitting given by 

X N = - 1 ~ T7 or TT ■ (3-12) 

m Nl MM y 1 

For instance, if M = 10 TeV and \il = ^ir — M 2 /Mx, one then finds xn w 10~ 12 - 
10 -11 . As we will see in Section 6, such small values of xn can lead to a resonant 
enhancement of CP asymmetries in the heavy Majorana neutrino decays. 

To obtain the sufficient and necessary conditions of CP invariance for any flavour 
structure of the one-generation model with two isosinglct neutrinos, one should 
use a more general approach, based on generalized CP transformations 47 for the 
fermionic fields: 

L L -> e^(L L f , v Rl -> VijivRjf, (3.13) 

where V is a 2 x 2 dimensional unitary matrix. Notice that the transformations 
given by Eq. (3.13) satisfy the SM symmetry of the mass- independent, conformal 
invariant part of the Lagrangian; only M v breaks this symmetry softly. In such an 
approach, 45,47 one looks for all possible weak-basis independent combinations that 
can be formed by Yukawa couplings and the neutrino mass matrix M v ', and are 
simultaneously invariant under the transformations (3.13). In this way, we find the 
condition 

lm r Tr{h)hM^M v M^h T h*M v ) = m Nl m N2 (m 2 Ni - m 2 N2 )\m{h n h* l2 f = 0, 

(3.14) 

where h = (tin , hi 2 ) is a row vector that contains the Higgs Yukawa couplings defined 
in the mass basis of isosinglet neutrinos. From Eq. (3.14), one readily observes that 
CP invariance holds if mjv x = mN 2 and/or one of the isosinglet neutrinos is massless. 
These considerations may be extended to models with til weak isodoublets and tir 
neutral isosinglets. In this case, there exist many conditions analogous to Eq. (3.14), 
which involve high-order terms in the Yukawa-coupling matrix h. However, not all 
of the conditions are sufficient and necessary for CP invariance. If we assume that 
Higgs triplets are not present in the theory, the total number of all non-trivial 
CP-violating phases is Ncp = ^l(«_r — l)- 48 
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Fig. 2. One-loop graphs contributing to the renormalization of the couplings \ INi, x° u l^i an d 
HviNi. 
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4. Renormalization 

At the tree level, CP violation in particle decays amounts to CPT violation and 
therefore vanishes identically. A non-vanishing contribution to CP asymmetries 
only arises at the one- loop level, considering the diagrams shown in Fig. 2. For this 
reason, it is important to discuss how one-loop renormalization applies to heavy 
Majorana-neutrino models 49 and its possible consequences on CP asymmetries. 

We start our discussion by expressing all bare quantities in terms of renormalized 
ones in the following way: 

"IL = E (<*"' + 2 5Z »T l ' L > l ° L = E (*"' + ^"'K' (41) 
v=i ' i'=i 

where unrenormalized kinematic parameters and fields are indicated by a super- 
script '0'. Note that 5Z$ collectively represents the wave- function renormalization 
constants of all components of the Higgs doublet $ (or $), i.e. the fields i X° an d 
iJ. In Appendix A, we give analytic expressions for Higgs and fermion self-energies. 
From these, one can easily see that the divergent part of the Higgs wave-function 
renormalization is exactly the same. In fact, 5Z& is universal in the limit Mh — > 0. 

Let us now consider that all quantities in the Lagrangian (3.7) are bare and we 
can substitute Eqs. (4.1) into that bare Lagrangian. In addition to the renormalized 
Lagrangian, which has the same structural form as the bare one, we then find the 
counter-term (CT) Lagrangian 

1 1=1 j=l V ni ' l' = l k=l ' 

(4.2) 

where L; = {vil, Il) T and SZ L = (SZ,SZ V ). Owing to charge and hyperchargc 
conservation on the vertices, it is not difficult to show by naive power-counting that 
the one-loop vertex corrections in Fig. 2(a)-(c) are ultra-violet (UV) finite (see also 
Appendix A). 

Despite the fact that vertex corrections are UV finite by themselves, the wave- 
function renormalizations of the Higgs and lepton isodoublets and that of neutrino 
isosinglets contain UV divergences that do not cancel. In accordance with the CT 
Lagrangian (4.2), one may require that all UV terms are to be absorbed into the 
definition of hij, i.e. 

1 til n R 

Sh l3 = - - ( h tj 6Zt + £ h Vj 8Zft + E hi k 6Z% ) . (4.3) 

i'=i fe=i 

It is important to stress that one-loop renormalization involves the dispersive parts 
of self-energies and effectively leads to a redefinition of the kinematic parameters, 
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whereas all absorptive corrections remain unaffected. Even though there might be 
some high-order dependence due to the choice of different renormalization schemes, 
we carry out the mass renormalization in the on-shell (OS) scheme. 50 As we will see 
in Section 7, this scheme has some field-theoretic advantages over other schemes, 
when applied to the resummation approach describing the mixing of two unstable 
particles. 

5. Resummation approach to unstable-particle mixing 

The consistent description of unstable particles within the conventional framework 
of perturbative S-matrix theory is an issue related to a number of field-theoretic dif- 
ficulties. Since unstable particles decay exponentially with time, they cannot appear 
as asymptotic in or out states in a process. Furthermore, finite-order perturbation 
theory breaks down. The usual propagator describing the unstable particle in the 
intermediate state of a given process displays a physical singularity when the parti- 
cle comes on its mass shell. One is therefore compelled to use resummation methods 
that treat unstable particles and unstable-particle mixing in a consistent way; this 
is a rather subtle issue within the context of gauge theories. 33 ' 30 

In a simple scalar theory with one unstable particle, Veltman 5 was able to 
show that, even if one removes the unstable particle from the initial and final states 
and substitutes it in terms of asymptotic states, the so-truncated S-matrix theory 
will still maintain the field-theoretic properties of unitarity and causality. Veltman's 
truncated S-matrix theory is rather useful to describe resonant processes in collider 
experiments where the initial and final states can be well prepared and detected. 
However, this formalism cannot directly be applied to the early Universe, as it does 
not take account of the many decoherentional collisions that an unstable particle 
may undergo with the thermal background before it decays. Therefore, one must 
seek a method that isolates the incoherent 2 part of an S-matrix amplitude. The 
new resummation method should include finite width effects in the mixing and de- 
cay of unstable particles. This will be done in an effective manner, by employing 
a procedure related to the Lehmann-Symanzik-Zimmermann formalism (LSZ). 53 
Then, the incoherent decay amplitude derived with this method may equivalcntly 
be embedded into a transition element 33 > 30 in line with Veltman's S-matrix formu- 
lation. As we will see in Section 9, the squared resummed decay amplitudes thus 
obtained will become the relevant collision terms entering the Boltzmann equations 
for the thermodynamic evolution of the Universe. 

We shall now demonstrate the effective resummation approach to unstable par- 
ticle mixing. Let us consider a theory with two neutral unstable scalars, e.g. S\ 
and S2 ■ The approach can then be extended to the case of unstable fermions such 
as heavy Majorana neutrinos. The bare (unrenormalizcd) fields Sf and their re- 
spective masses M° may then be expressed in terms of renormalized fields Si and 
masses Mi in the following way: 

Sf = zlj 2 S 3 = (Sij + \5Z ij )s j , (5.1) 
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lim 




?- 1/2T A"V) 



Fig. 3. Diagrammatic representation of the renormalized n — 1-non-amputatcd amplitude, Si,..., 
and the LSZ reduction formalism. 



(Mf) — Mf + 5Mf 



(5.2) 



Here and henceforth, summation is understood over repeated indices that do not 

1 /2 

appear on both sides of an equation. In Eqs. (5.1) and (5.2), ZJ- and SMi are 
the wave-function and mass renormalization constants, respectively, which can be 
determined from renormalization conditions imposed on the two-point correlation 
functions, ILj(p 2 ), for the transitions Sj — > Si in some physical scheme, such as 
the on-mass-shell (OS) renormalization scheme. 50 More details may be found in the 
appendix. 

In order to include the mixing of the unstable scalars, we must first calculate 
all the SiSj Green functions, with i,j = 1, 2. After summing up a geometric series 
of the self-energies Hij(p 2 ), the full propagators may be obtained by inverting the 
following inverse propagator matrix: 



A-V) 



p 2 - (Mf) 2 + imp 2 ) 
iW) 



ni 2 (p 2 ) 

(M 2 ) 2 + Ii 22 {p 2 ) 



The result of inverting the matrix in Eq. (5.3) may be given by 

n 2 2 ( P 2 ) 



Ah(p 2 ) 
A 22 (p 2 ) 
Ai 2 (p 2 ) 



p 2 (Af 1 °) 2 + n 11 ( P 2 )- 
p 2 - (M 2 °) 2 + n 22 (p 2 )- 
a 21 ( p 2 ) = -n 12 ( s ) 



P 2 -(M 2 °) 2 + n 22 (p 2 ) 
n 2 2 ( P 2 ) 



p2_ (M )2+nii(p2) 

(p 2 -(M 2 °) 2 + n 22 ( P 2 )) 



(5.3) 

(5.4) 
(5.5) 



^-(M^ + iiMp 2 )) -tiUp 2 ) 



(5.6) 



where ITi 2 (p 2 ) = II 2 i(p 2 ). Moreover, we observe the crucial factorization property 
for the off-diagonal (i ^ j) resummed scalar propagators 

iW) 



A^V) 



A u (p 2 ) 



p 2 - (m 7 ) 2 + n w ( P 2 ) 
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n tJ ( P 2 ) 

P 2 - (M°) 2 + U u (p 



2_r^2^TT..^ A «V) • (5-7) 



The resummed unrcnormalizcd scalar propagators (p 2 ) are related to the respec- 
tive renormalized ones Ay(p 2 ) through the expression 

A«(P 2 ) = ^ 2 A m „(p 2 )^f T , (5.8) 

where Aij(p 2 ) may be obtained from Eqs. (5.4)-(5.6), just by replacing Mf with 
Mj and Hy(p 2 ) with n i:) (p 2 ). Note that the property given in Eq. (5.7) will also 
hold true for the renormalized scalar propagators Ay(p 2 ). 

Suppose that we wish to find the effective resummed decay amplitude 7^ for 
the decay Si to n light stable scalars Si t , . . . , Si n . In analogy to the LSZ formalism, 
one starts with the Green function describing the transition shown in Fig. 3, and 
amputates the external legs by their inverse propagators. For the stable external 
lines Si t , . . . , Si n , the procedure is essentially the same with the usual LSZ formal- 
ism. This formalism may then be extended to the external line describing the SiSj 
system. The intermediate steps of this procedure are given by 

P i 

= lim 

p 2 -*M? 



rj-amp yl/2 /j-amp 7 l/2 Hmt {P )(1 ^mi) 

p 2 -mi + n mro ( P 2 ) 

%(M 2 )(l-%) 
M 2 -M 2 + n w (M 2 ) 

where and Tj... are the renormalized transition elements evaluated in the stable- 
particle approximation. One should bear in mind that the OS renormalized self- 
energies Tlji(Mf) in Eq. (5.9) have no vanishing absorptive parts, as renormalization 
can only modify the dispersive (real) part of these self-energies. The reason is 
that the CT Lagrangian must be Hermitian as opposed to the absorptive parts 
which are anti-Hcrmitian. In fact, these additional width mixing effects are the 
ones we wish to include in our formalism for decay amplitudes and are absent 
in the conventional perturbation theory. It is also important to observe that our 
approach to decays is not singular, i.e. Si... displays an analytic behaviour in the 
degenerate limit M 2 — > M 2 , because of the appearance of the imaginary term 
ilmJIjj(Mf) in the denominator of the mixing factor present in the last equality 
of Eq. (5.9). Finally, we must stress that the inclusion of these phenomena has 
been performed in an effective manner. Since the decaying unstable particle cannot 
appear in the initial state, 51 the resummed decay amplitude must be regarded as 
being a part which can effectively be embedded into a resummed S-matrix element. 33 
This resummed S-matrix element describes the dynamics of the very same unstable 
particle, which is produced by some asymptotic states, resides in the intermediate 
state, and subsequently decays either directly or indirectly, through mixing, into 
the observed final states. 
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The resummation approach outlined above can now carry over to the mixing 
between two unstable fermions, call them /i and f 2 . As we did for the case of 
scalars, we express the bare left- and right-handed chiral fields, and f Ri (with 
i = 1, 2), in terms of rcnormalizcd fields as follows: 

f° 

J Li 



fm — Z Rij fRj i 



(5.10) 



where Z^f? {Z R ?) is the wave- function renormalization constant for the left (right)- 
handed chiral fields, which may be determined from the fcrmionic self-energy tran- 
sitions fj — > fi, Ejj(j/), e.g. in the OS renormalization scheme. 49 Analogously to 
Eq. (5.3), the resummcd fcrmion propagator matrix may be obtained from 

-l 



(5.11) 



where m\ 2 are the bare fermion masses, which can be decomposed into the OS 
renormalized masses mi l2 and the CT mass terms £7711,2 as 2 = mi l2 + <5mi, 2 . 
Inverting the matrix-valued 2x2 matrix in Eq. (5.11) yields 

1 i-i 



S11OQ 

Sl2(j0 

5 2 i (jO 



m° 2 + E 22 (^) 



Si 2 (jO 



/f- m% + E 22 ( 
1 



Sn(lO 



£ 2 i(0 

Ei2(jQ 



= -5 u (^£i2(^) J* - m% + S 22 (^) 



£l2 (?0 5*22 (/0 , 

S 2 2(jOS 2 i(iO ™i + £iiW] _1 
* - m° 2 + £ 22 (jQ1 _1 £21 OOSiikO ■ 



,(5.12) 
L ,(5.13) 



(5.14) 



(5.15) 



Equations (5.14) and (5.15) show that the resummed propagators S\ 2 (j?) and S 2 \(jf) 
are endowed with a factorization property analogous to Eq. (5.7). Similarly, the 
renormalized and unrcnormalizcd resummcd propagators are related by 



2 Rim P r) ^« 



>(?9 (z_ 



l/2f 
Lnj 



Pe 



Zr'^Pl) 



(5.16) 



where the caret on SV, (jf) indicates that the resummed fermionic propagators have 
been renormalized in the OS scheme. Moreover, the renormalized propagators Sij 
(jf) may be obtained by Sij(ff) in Eqs. (5.12)-(5.15), if the obvious replacements 
m® — > rrii and Ejj(^) — > Ey(^) are made. 

By analogy, one can derive the resummed decay amplitude, of the unstable 
fermion — > A, as we did for the scalar case. More explicitly, we have 

,1/2 



"Ti... 



1 k... \ L . 



1/2 
Lkm 



Pl + z R <: m p r ) s n M (zffi p r + Pl) 



x(ZI^Pr 



^•/ 2t P i )s« 1 (iO«i(p) 



= 3. 



(5.17) 



(1 - Si^Tj... Ej 



m. 



+ £jj(/0 
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Fig. 4. e- and s'-type CP violation in the decays of heavy Majorana neutrinos. 

Again, 7$... represent the respective renormalized transition amplitudes evaluated in 
the stable-particle approximation. The amplitudes %... also include all high-order 
n-point functions, such as vertex corrections. Based on the formula (5.17), we shall 
calculate the CP asymmetries in the decays of heavy Majorana neutrinos in the 
next section. 



6. CP asymmetries 

The resummation approach presented in the previous section may be applied to 
describe e- and e'-type CP violation in heavy Majorana neutrino decays shown in 
Fig. 4. The same formalism may also be used to determine the collision terms for 
the inverse decays, which occur in the formulation of the Boltzmann equations (see 
also Section 8). 

Let us consider the decay N\ — > l~x + m a model with two right-handed neutri- 
nos. The inclusion of all other decay channels is then obvious. We shall first write 
down the transition amplitude responsible for e-type CP violation, denoted as T^f 1 , 
and then take CP-violating vertex corrections into account. Applying (5.17) to 
heavy Majorana neutrino decays, we obtain 

-l 



z 



(e) 



fin uiP R u Nl - iha uiP R i> - m N2 + iZ%2 S (f>) 



where the absorptive part of the one-loop transitions Nj 
the general form 



ESf-kQuiv,, (6.1) 
Ni, with i,j = 1, 2, has 



with 



MAv 2 ) = 



hi'iht, 



3 1 
2 + 2 



p2 



(6.2) 



(6.3) 



32tt 

In the limit M H -» 0, Eq. (6.3) gives A tJ = h lH hp } / (16tt) . The CP-transform 
resummed amplitude describing the decay N\ — » l + x~ , T~n, , reads 



1 N, 



= h* a v Nl P L vi 
= h* n uiP L u Nl 



ih l2 v Nl Z? 2 s (- f) 
ih* l2 uiP L f - m N2 



-l 



iT ahs ( 



—abs 



Plvi 
(6.4) 
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where 

S7(jf) = Aijtf) i»P R + A*^) $P L (6.5) 

is the charge-conjugate absorptive self-energy. The last step of Eq. (6.4) is derived 
by making use of the identities 

u(p,s) = Cv T (p,s), C^C'- 1 = - 7 J. (6.6) 

The expressions in Eqs. (6.1) and (6.4) may be simplified even further, if the Dirac 
equation of motion is employed for the external spinors. Then, the two resummed 
decay amplitudes, and , take the simple form 

-.(e) - D T, m%(l + iA 22 )A* 21 +m Nl m N2 A 21 1 

1 - u l P R u Nl [h ll - l h l2 m%i(1 + lA22)2 _ mk J, (6-7) 



T^l = uiP L u Nl h* n - ih 



H2 ' 



m 2 Ni (l + iA 22 )A 21 + m Nl m N2 A* 21 



m 2 Ni (l +iA 2 2) 2 - m 



N 2 



(6.8) 



In addition, the respective transition amplitudes involving the decays N 2 — > I \ + , 
T^ 2 \ and iV 2 — ► l + x ~ ■> ^n 2 -> ma y obtained by interchanging the indices '1' and 
'2' everywhere in Eqs. (6.7) and (6.8). 

In order to study the e- and e'-type mechanisms of CP violation in heavy Ma- 
jorana neutrino decays, we define the following CP-violating quantities: 



e Ni = — f- _ ' , fori = 1,2, (6.9) 

i4?i 2 + ra 2 

e N = 1n ; +1n ; — ^ — ^ff. (6.io) 
l^l 2 + l^l 2 + TO 2 + ra 2 

Correspondingly, the CP-violating parameters e' N . and e' N may be defined by 



iri r £,) i 2 - ir^ ' 2 



ek = M^TT =§r> for z = 1,2, (6.11) 

\T^\ 2 + \T% ] \ 2 



'N 



(e')|2 _l |-r(e')|2 _ , _ . 

(6.12) 



i^ 1 ) i 2 + i^: , i 2 -i^ 1 t-i^ 2 



i< 5 i 2 + i4:°i 2 + iO + 1^?! 2 



The last parameters quantify CP violation coming exclusively from the one-loop 
irreducible vertices. In Eqs. (6.9) and (6.10), the parameters ejv 4 and £n share 
the common property that they do not depend on the final state that N decays, 
despite the fact that the individual squared matrix elements do. In general, both 
s- and e'-type contributions are not directly distinguishable in the decay widths 
T(Ni — > ^x*), unless » e' N . and vice versa, for some range of the kinematic 
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parameters. Evidently, the physical CP asymmetries are given by 



fori = 1,2, (6.13) 



where L refers to all fcrmionic degrees of freedom of the leptonic isodoublet that 
heavy Majorana neutrinos can decay. Nevertheless, the parameters e^., en, £jy. 
and e' n defined above are very useful to determine the contributions due to the 
different mechanisms of CP violation. 

We now turn to the calculation of the CP-violating contribution, which is entirely 
due to the heavy-neutrino self-energy effects. Substituting Eqs. (6.7) and (6.8) into 
(6.9), we arrive at the simple formulas 21 

_ ImjJ^hnf Am 2 N m Nl T N2 



Im(fe ; * 1 fe; 2 ) 2 Am 2 N m N2 r Nl 
V \h n \*\h a \* (Am 2 N ) 2 + m%T 2 Ni ' 



where Am^, = m 2 Ni — m% 2 and 

T Ni = ^ m Nt (6.17) 

are the decay widths of the heavy Majorana neutrinos. Equations (6.15) and (6.16) 
are a very good approximation for any range of heavy- neutrino masses of interest. 
Both CP asymmetries £at x and en 2 are of the same sign and go individually to zero 
when Amjsj — > 0, as it should be on account of Eq. (3.14). In the conventional 
perturbation theory, the width terms rn 2 N Y 2 N2 and m 2 N2 Y 2 Ni occurring in the last 
denominators on the RHS of Eqs. (6.15) and (6.16) are absent. This very last fact 
is precisely what causes a singular behaviour when the degeneracy between the two 
heavy Majorana neutrinos is exact. On physical grounds, however, the only natural 
parameter that can regulate such a singularity is the finite width of the heavy 
neutrinos, which is naturally implemented within the resummation approach. 

From Eqs. (6.15) and (6.16), it is not difficult to derive the sufficient and nec- 
essary conditions for resonant enhancement of CP violation. To be specific, CP 
violation can be of order unity if and only if 

r r 

(i) m Nl - m N2 ~ ±A 22 m N2 = and/or A 11 m Nl = , (6.18) 

llm^M 2 ! -i , fi1Q , 

Scp = M 2 IM 2 " 1 • (6 - 19) 

Before we present numerical estimates of CP asymmetries, we calculate for 
completeness the contributions to CP violation arising entirely from vertex effects. 
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The e'-type contributions can be significant for large differences of heavy neutrino 
masses, e.g. for mjvi — mN 2 <~ or itin 2 . In this regime, both e-type and e'-type 
effects are of comparable order. 19 It is useful to define first the function 



F(x,a) = - a - (l + ^ln^ 1 " + ^ ) 



(6.20) 



With a — 0, F(x, a) reduces to the Fukugita-Yanagida loop function f(x) = y/x[l — 
(1 + .x)ln(l + 1/x)]. 15 Then, L-violating absorptive parts of the one-loop vertices 
X + lNi, x°viNi and HviNi, shown in Figs. 2(a)-(c), are given by 



32tt Jp 1 



i>Ph f 



( Nj 
V p 2 







(6.21) 



i>p L 



Here, we have assumed that the external decaying heavy Majorana neutrinos are 
off-shell, whereas the leptons and Higgs fields are on their mass shell. The complete 
analytic expressions are calculated in the appendix. Using Eqs. (6.21) and (6.22) 
and neglecting wave-function contributions, we compute the e'-type CP asymmetry 
in the conventional perturbation theory. Considering all decay channels for the 
decaying heavy Majorana neutrino, e.g. N\, we find 



Im^M 2 



l^\hii\ 2 [l + k{l-M%/m? Ni Y] 



f 5 p( m N 2 



N% N% "\2Vi 



m 



£)]}• 



Ill- 



No 



M 



III: 



N 1 



Til 



(6.23) 



In the limit Mr — > 0, the last formula simplifies to the known result 15,16,17,18 



Nl 8tt|^i| 2 



(6.24) 



Unlike ejvi, c'n does not vanish in the degenerate limit of the two heavy Majorana 
neutrinos N\ and A^. However, when the value of approaches that of m.jv 2 , the 
e'-type part of the transition amplitude squared for the N\ decay becomes equal but 
opposite in sign to the respective one of the N 2 decay. As a result, these two e'-type 
terms cancel one another, leading to the vanishing of the CP-violating parameter 
e' N defined in Eq. (6.12). Consequently, as opposed to e effects, e' effects cannot 
become resonant for any kinematic region of mass parameters. 

Both £ and e' contributions can be included into the resummed decay ampli- 
tudes. Considering Eqs. (6.21) and (6.22) into account, we obtain 



Tr 



UlPl 



[hi 



H2 
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(6.25) 



x i> - m N2 + iS 22 [ft) H 21 {f)\u Nl , 



(6.26) 



where the notation of the off-shell one- loop vertices has been simplified to Vjf s (jf). 
The vertex functions V^* '(ft) are the charge conjugates of Vf^ s { f) and may hence 
be recovered from Eqs. (6.21) and (6.22), by taking the complex conjugate for the 
Yukawa couplings and replacing Pp with Pj, . Although the calculation of the CP- 
violating observables 8n { defined in Eq. (6.13) is quite straightforward from Eqs. 
(6.25) and (6.26), it is not very easy to present analytic expressions in a compact 
form. 

To gauge better the dependence of CP asymmetries on the heavy neutrino 
masses, we shall adopt two simple scenarios with two-right handed neutrinos that 
mix actively with one lepton family I only: 



We assume that N 2 is always heavier than N\, i.e. mjvi < mjv 2 - The above two sce- 
narios comply qualitatively with Sakharov's third requirement of out-of-equilibrium 
condition (see also discussion in Section 8). In view of Eq. (6.19), both scenarios I 
and II given above represent maximal cases of CP violation with 5c p = 1- There- 
fore, results for any other model may readily be read off by multiplying the CP 
asymmetries with the appropriate model-dependent factor Sep- 

Figure 5 exhibits the dependence of the CP asymmetries as a function of the 
parameter xn for scenario I. The parameter xn defined in Eq. (3.12) is a measure 
of mass degeneracy for the two heavy Majorana neutrinos Ni and N 2 - We divide 
the range of values for the parameter xn into two regions: the first region is plotted 
in Fig. 5(a) and pertains to the kinematic domain where resonant CP violation due 
to heavy-neutrino mixing occurs. The second one, shown in Fig. 5(b), represents 
the kinematic range far away from the resonant CP-violating phenomenon. The 
dotted line in Fig. 5(a) gives the prediction of en 17 when Eq. (6.15) is calculated 
in the conventional finite-order perturbation theory. Obviously, e p ^ diverges for 
sufficiently small values of xn, e.g. xm < 1CP 13 . If resummation of the relevant 
fermionic self-energy graphs is considered, the prediction for £jvi becomes analytic 
and is given by the dashed line in Fig. 5. The Sn 1 line shows a maximum for 
xn ~ 10~ 13 . In agreement with the conditions in Eqs. (6.18) and (6.19), CP 
violation may resonantly increase up to order unity. 36,30 The solid line in Fig. 5 
displays the dependence of the CP-violating parameter Sn in Eq. (6.14) on xn, 
where e'-type contributions are included. The latter are very small in this scenario, 
so as to account for the BAU, e.g. e' N ~ 1CP 16 . Finally, we comment on the fact 
that 8n vanishes in the CP-invariant limit x n — ► 0, as it should be on account of 



I. m Nl = 10 TeV, 
II. m Nl = 10 9 TeV, 



h n - 10- 
h n = 10 



6 



h n = 10~ 6 (l + i), 
h n = 10~ 2 (^). (6.27) 



-2 



Eq. (3.14). 
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Figures 6 and 7 give numerical estimates of CP asymmetries in scenario II. The 
difference of this model with that of scenario I is that the e'-type effects may not 
be negligible in the off-resonant region, as can be seen from Figs. 6(a) and 7. In 
particular, for values of the parameter xn < 10~ n or xn > 1, the individual e' N - 
and £^- 2 -type contributions may prevail over the £-type ones. Models with xn > I 
have been extensively discussed in the literature. 15,16,17,18 Numerical estimates for 
such models are displayed in Fig. 7. We first focus our attention on the domain 
with xn < 1CP 2 . In Fig. 6(a), we observe that e' N and e' N , represented by the 
dotted lines, do not vanish in the CP-invariant limit ij/ — > 0, as opposed to e^. 
As a consequence, the CP asymmetry Sn 1 in Eq. (6.13), in which both e^- and 




Fig. 5. Numerical estimates of CP asymmetries in scenario I. 
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Fig. 6. Numerical estimates of CP asymmetries versus mN 2 / m N 1 — 1 in scenario II. 



£ jvrtyP e terms are considered within our formalism, docs not vanish either. The 
reason is that the physical CP-violating parameter in this highly degenerate mass 
regime for Ni and N 2 is the observable Sn defined in Eq. (6.14). In fact, Sn and en ± 
share the common feature that both tend consistently to zero as x n — > 0. This fact 
must be considered to be one of the successes of the resummation approach. Again, 
CP violation is resonantly amplified, when the condition in Eq. (6.18) is satisfied, 
as can be seen from Fig. 6(b). Finally, we must remark that — Sn flips sign and 
eventually becomes negative for xn ^> 1, as can be seen from Fig. 7. However, in 
this kinematic range, we must consider a further refinement into the definition of 
Sn- The effect of the different dissipative Boltzmann factors multiplying the decay 



Heavy Majorana neutrinos and baryogenesis 25 

rates of the heavy Majorana neutrinos N\ and N 2 must also be included in 5jy. 
These phenomena will be taken into account in Section 8. 

7. Unitarity and CPT invariance in the resummation approach 

It is interesting to see how the resummation approach preserves CPT invariance 
and unitarity. 30 An immediate consequence of unitarity and CPT symmetry is that 
CP violation in the L-violating scattering process — > L c $ is absent to order 
^6 54,55 w jjj concentrate on the resonant part of the amplitude, as it is the 
dominant one. 

Our aim is to show that to 'one loop', 

Acp = / dLIPS |T L r |U L c $ | 2 - J dLIPS IT^^tl 2 = 0, (7.1) 

where LIPS stands for the two-body Lorentz-invariant phase space. For simplic- 
ity, we omit external spinors and absorb irrelevant constants in the definition of 




Fig. 7. Numerical estimates of CP asymmetries as a function of mjvj/m^ — 1 in scenario II. 
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the Yukawa-coupling matrix h = (hn,hi2). Using matrix notation, the resummed 
transition amplitudes are written 

^ s »^ c$ = hP R S(^P R h T , TTc^ l ^ = h*P L S{f)P L h), (7.2) 

with S(jf) — S T {jf) being the CP/T-conjugate propagator matrix of S(jf). In writing 
the CP/T-conjugate amplitude T jcj^^^t , we have employed the identities (6.6) 
for spinor objects and made use of the rotational symmetry of the amplitude. The 
latter has the effect of reversing the spatial components of the four momenta. We 
also neglect possible P-odd spin correlations involving external leptons since they 
will be averaged away when forming the matrix element squared. 

We start the proof by noticing that as a consequence of CPT invariance, 

ITt*t_>.i,*t I = I^l°*->.e,c*I ■ (7.3) 

This equality is indeed valid, since 

\hP R SMP L rf\ = \h*P L S T (tf)P R h T \ = \h*P L S(^)P R h T \ . (7.4) 

Unitarity of the theory prescribes the following relation governing the resummed 
propagators: 

S^ifi) - S~^(f) = -i J rfLIPS i>(h T h*P L + h)hP R ) . (7.5) 

This last relation is also known as the optical theorem. Based on the optical theo- 
rem, we can prove the equality 

J dLIPS \Tl^^ L ^ L c^\ 2 = J cJLIPS \T L c^ L ^ tL c^\ 2 . (7.6) 

Indeed, using Eq. (7.5), we find 

= j dLIPS hP R S(f) i>{h T h*P L + h)hP R )SH$)P L h) 
= -ihP R [S{f) - &{fj\P L rf = 2hP R lmS{tf)P L rf , (7.7) 
and for the CP-conjugate total rate, 

J dUPS\Tl c c %^ L<s>f L c*\ 2 = 2h*P L lmS(fi)P R h T 

= 2hP R lmS(^)P L hK (7.8) 

As the RHSs of Eqs. (7.7) and (7.8) are equal, the equality (7.6) is obvious. Sub- 
tracting Eq. (7.3) from Eq. (7.6), it is not difficult to show that Acp vanishes at 
the one-loop resummed level. We should remark that the resummation approach 30 
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satisfies CPT and unitarity exactly, without recourse to any re-expansion of the 
resummed propagator. If we also include resummed amplitudes subleading in the 
Yukawa couplings, then residual CP-violating terms that are formally of order 
and higher occur in Aqp- These terms result from the interference of two resummed 
amplitudes containing one-loop vertex graphs. Because of unitarity, however, the 
residual CP-violating terms of order hf t and h}® will cancel at two loops together 
with respective CP-violating terms coming from one-loop 2 — > 4 scatterings, and so 
on. 

In the approach under consideration, 30 the physical transition amplitude is ob- 
tained by sandwiching the resummed propagators between matrix elements related 
to initial and final states of the resonant process. Therefore, diagonalization of 
S(]i) is no longer necessary, thereby avoiding possible singularities emanating from 
non-diagonalizable (Jordan-like) effective Hamiltonians [or equivalently S^ 1 (]/)]. 30 
In fact, such effective Hamiltonians represent situations in which the CP-violating 
mixing between the two unstable particles reaches its maximum and physical CP 
asymmetries are therefore large. In such a case, the complex mass eigenvalues of 
the effective Hamiltonian are exactly equal. 

To see this point in more detail, let us consider the following effective Hamilto- 
nian for the NiN 2 system: 

rriN + a — i\b\ —ib 

—ib* raw — a — i\b\ ' 

(7.9) 

in the approximation ft — > w mi w m 2 . In Eq. (7.9), the parameters a 
and b are real and complex, respectively, and mi = + a, m 2 = m^ — a. 
The complex parameter b represents the absorptive part of the one-loop neutrino 
transitions Ni — > Nj . Unitarity requires that the determinant of the absorptive part 
of H(tf) be non-negative. For the effective Hamiltonian (7.9), the corresponding 
determinant is zero. One-generation models naturally lead to such an absorptive 
effective Hamiltonian. If a = \b\, the two complex mass eigenvalues of Tt($) are 
exactly degenerate and equal to m^v — i\b\. Then, the effective Hamiltonian cannot 
be diagonalizcd via a similarity transformation in this limit, i.e. the respective 
diagonalization matrices become singular. 30 

An interesting question one may raise in this context is the following. Since 
models with non-diagonalizable effective Hamiltonians lead to an exact equality 
between their complex mass eigenvalues, how then can this fact be reconciled with 
the CP-invariance condition (3.14)? According to the condition (3.14), any effect 
of CP violation must vanish identically, and should not even be large! To resolve 
this paradox, one should notice that in the presence of a large particle mixing, the 
mass eigenstates of S ,_1 (y) are generally non-unitary among themselves, whereas 
the OS-renormalized mass eigenstates 49 form a well-defined unitary basis (or any 
other renormalization scheme that preserves orthonormality of the Hilbert space), 
upon which perturbation theory can be formulated order by order. Therefore, the 
field-theoretic OS renormalizcd masses are those that enter the condition of CP 
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invariance given by Eq. (3.14). Consequently, if the two complex mass eigenvalues 
of the effective Hamiltonian are equal, this does not necessarily entail an equality 
between their respective OS renormalized masses, and therefore absence of CP 
violation as well. 30 

8. Boltzmann equations 

The thermodynamic evolution of the system in the radiation-dominated era of the 
Universe may be described by a set of coupled Boltzmann equations (BE's). 4 ' 5,2 ' 38 
These equations determine the time evolution of the lepton-number asymmetry 
which will be converted into the observed BAU by sphalerons. We shall solve the 
BE's numerically and present results for the expected BAU within the two different 
democratic-type scenarios I and II discussed in Section 6. Finally, we will give 
estimates of the finite-temperature effects, and discuss their impact on resonant CP 
violation through mixing. 

Before solving numerically the BE's, it is instructive to discuss first the out- 
of equilibrium constraints on heavy neutrino decays. Sakharov's third necessary 
condition requires that the expansion rate of the Universe be smaller than the 
decay rate of any L-violating process. The most dominant L-violating process is 
the decay of the heavy Majorana neutrinos themselves, (cf. Eq. (6.17)). To a 
good approximation, we have the approximate inequality 

r Ni (T = m Ni ) < 2KH(T = m Ni ), (8.1) 

where K w 1-1000 is a factor quantifying the deviation of the decay rates from the 
expansion rate of the Universe, and H{T) is the Hubble parameter 

H(T) = 1.73 gl' 2 ^—- ■ (8.2) 

Mpi anc k 

with Mpianck = 1.2 x 10 19 GeV and g„ w 100-400 being the number of active degrees 
of freedom in usual extensions of the SM. Then, the out-of equilibrium constraint 
in Eq. (8.1) translates into the bound 

|/^| 2 £ 7.2XxlO- 14 ( I ^). (8.3) 

Although not mandatory, this very last constraint may be applied to all Yukawa 
couplings. 

As we have discussed in Section 2 (see also Eq. (2.13)), above the electroweak 
phase transition the B + L-sphalcron interactions convert approximately one-third 
of the lepton-to-entropy density ratio Yl = til/s into a baryon-to-entropy density 
ratio Y B = n B /s, i.e. 11,38 

The last approximate equality represents the asymptotic solution of the relevant 
BE's. 16 ' 21 From Eq. (8.4), we see that Yg can be in the observed ball park, i.e. 
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Yb ~ 1CP 10 , if \dNi\/K are of order 10 _7 -10 -6 . Clearly, CP asymmetries of order 
unity allow for very large values of K. As a consequence, the thermal plasma 
can then be rather dense and the conditions of kinetic equilibrium in BE's can 
comfortably be satisfied even within the minimal leptogenesis scenario under study. 

We now turn to the discussion of BE's. The lepton asymmetry for a system with 
two heavy Majorana neutrinos is determined by the coupled system of BE's 2 ' 16 



h dHjlN, = - - 1 



l Ni 



+ 3Hn Ni = -(J± -l) lNi (8.5) 



dn L , „„ \ x (n Ni \ n L 1 n L 

~w + 3HnL = ^ I** Uf " l ) ~ 2^\ 1N > ~ (8 - 6) 

where njVj , nL = ni — nj are the densities of the number of Ni and the lepton- number 
asymmetry, respectively, and n e £. and n e ] q are their values in thermal equilibrium. 
The Hubble parameter H = (dR/dt)/R determines the expansion rate of the Uni- 
verse and also depends on the temperature T, through the relation in Eq. (8.2). 
In Eqs. (8.5) and (8.6), r jN i and -f a are the decay and scattering collision terms, 
respectively: 

IN, = ^ i K2{m%JT) T Ni , (8.7) 



(8.E 



rp p OO 

la = — J dss^K^s/T)a'{s). 
Here, sthr is the threshold of a generic process a + b — > c + d, and 

ct'(s) = \ 9{^-m a -m b ) A(l, ^, ^) a(s) (8.9) 

with A(x, y, z) = (x — y — z) 2 — Ayz. In Eq. (8.7), K\{z) and K 2 (z) are the modi- 
fied Bessel functions defined in Ref. 56 The cross section <r(s) mainly comprises the 
scatterings L c $ — > and its CP-conjugate process — > and is eval- 

uated at T = by subtracting all those real intermediate contributions that have 
already been taken into account in the direct and inverse decays of heavy Majorana 
neutrinos. 5 The collision term 7 CT acts as a CP-conserving depletion term, which is 
formally of order 7^. at T k, mjv 4 - 2 There is also the AL = 2 reaction <-> LL, 
which is much weaker than the latter, as long as the out-of-equilibrium constraint 
on the Yukawa couplings in Eq. (8.3) is imposed. Finally, there exist additional con- 
tributions to the BE's, 16 coming from processes such as NiL <-» Qiia, NiQi <-> Ur. 
These contributions are quite strong at very high temperatures, T 3> Tnjv i5 and 
lead to a decoherence phenomenon between the heavy neutrinos N\ and N 2 . At 
the crucial leptogenesis epoch, when T m m^, the rates of the latter processes are 
kinematically suppressed and smaller than the decay rates of the heavy Majorana 
neutrinos. 2 
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Many applicable assumptions are involved in BE's (8.5) and (8.6). More de- 
tails may be found in Ref. 2 First, we have considered the Friedmann-Robertson- 
Walker model in the non-relativistic limit. Second, we have adopted the Maxwell- 
Boltzmann statistics, which is a good approximation in the absence of effects that 
originate from Bose condensates or arise from a degeneracy of many Fermi degrees 
of freedom. Third, we have assumed that the lepton and Higgs weak isodoublets, L 
and $, are practically in thermal equilibrium, and neglected high orders in n^ln\ q 
and (J/Vj- In this context, it has also been assumed that the different particle species 
are in kinetic equilibrium, i.e. that the particles may rapidly change their kinetic 
energy through elastic scatterings but the processes responsible for a change of the 
number of particles are out of equilibrium. These out-of-equilibrium reactions are 
described by the BE's (8.5) and (8.6). 

To solve these BE's numerically, it proves useful to make the following change 
of variables: 



m Nl 



t = 



1 



(8.10) 



T ' 2H(T) 2H(x = 1) ' 

Such an ansatz is also valid for the radiation-dominated phase of the Universe while 
baryogenesis takes place. Then, we define the parameters 



K 



Kijx) 

K 2 (x) H(x = l) 



7 = 



K^K^x) T N2 



Ki(x)K 2 (£x) T Nl ' 



^ , (8.H) 



with £ = mN 2 /rriN 1 > 1. In addition, we introduce the quantities Yjv 4 — n^Js and 
Yl = til/s, where s is the entropy density. In an isentropically expanded Universe, 
the entropy density has the time dependence s(t) — const, x R~ 3 (t) and may be 
related to the number density of photons, where g* is given after 

Eq. (8.2). Emplyoing the above definitions and relations among the parameters, we 
obtain the BE's in terms of the new quantities Y^, Yjv 2 and Yj,: 



dY Nl 

dx 
dY N2 

dx 
dY L 

dx 



= -(Y Nl -Y%)Ka?, 
- -(Y^-Y^Kx 2 



(8.12) 
(8.13) 



(Y Nl -Y%)5 Nl + (Y N2 -Y^)jS N2 - -g*Y L {Y e N \ +1 Y e N l) 



1N 1 



Kx 2 



(8.14) 



The heavy-neutrino number-to-entropy densities in equilibrium Y^.(x) are given 
by 5 

YnS x ) = T~ / dzzVz 2 -x 2 e- z = x 2 K 2 (x) , (8.15) 
°9* Jx °9* 

and Y^ q 2 (x) = Yj^(£x). The differential equations (8.12)-(8.14) are solved numeri- 
cally, using the initial conditions 



Yv,(0) = Yv 2 (0) = Y^(0) 



Y£(0) and Y L (0) = 0. 



(8.16) 
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These initial conditions merely reflect the fact that our Universe starts evolving from 
a lepton-symmetric state, in which the heavy Majorana neutrinos are originally in 
thermal equilibrium. Here, we should remark that the low-temperature limit of the 
numerical predictions does not strongly depend on the initial conditions (8.16), if L- 
violating interactions are in thermal equilibrium at T > m^r. . The reason is that at 
very high temperatures, the BE's (8.12)~(8.14) exhibit a running independent fixed 
point, and any initial value for Y/Vi, Yjy 2 and Yl is rapidly driven to the thermal- 
equilibrium values given by Eq. (8. 16). 57 After the evolution of the Universe to 
temperatures much below mjv i; a net lepton asymmetry has been created. This 
lepton asymmetry will then be converted into the BAU via the sphalerons. During 
a first order electroweak phase transition, the produced excess in L is also encoded 
as an excess in B, which is given by Eq. (2. 13). 11,38 The observed BAU is Y§ bs = 
(0.6 — 1) x 10~ 10 , 2 which corresponds to an excess of leptons ~Y£ bs w 10~ 9 — 1CP 10 . 
In the latter estimate, we have included the possibility of generating the BAU via 
an individual lepton asymmetry. 39 

Figure 8 shows the dependence of Yl{x) on x = m^ 1 jT for two representative 
scenarios defined in Eq. (6.27) for different values of xm — mN 2 /niN 1 — l. 21 The 
observed range for Yl, Y£ bs , is indicated with two confining horizontal dotted lines. 
In scenario I (Fig. 8(a)), a heavy-neutrino mass splitting xn of order 10 _9 -10 -8 
is sufficient to account for the BAU. For comparison, it is worth mentioning that 
the degree of mass degeneracy between Kl and K$ is of order 10~ 15 , which is by 
far smaller than the one considered here. We find that the e-type CP violation is 
dominant, whereas e'-type effect are extremely suppressed. Numerical estimates 
for the second scenario are displayed in Fig. 8(b). This scenario is closer to the 
traditional one considered in Ref. 15 Here, it is not necessary to have a high degree 
of degeneracy for N\ and N 2 to get sufficient CP violation for the BAU. In this case, 
both s- and e'-type mechanisms of CP violation are equally important. Therefore, 
the main consequence of £-type CP violation is that the leptogenesis scale may be 
as low as 1 TeV, even for models with universal Yukawa couplings. 21 

In the scenario of leptogenesis induced by mixing of heavy Majorana neutrinos, 
one may have to worry about effects, which could affect the resonant condition of 
CP violation in Eq. (6.18). For instance, there may be broadening effects at high 
temperatures due to collisions among particles. Such effects will contribute terms 
of order to the Ni widths and are small in general. 5,58 On the other hand, finite 
temperature effects on the T = masses of particles may be significant. Because 
of the SM gauge interactions, the leptons and the Higgs fields receive appreciable 
thermal masses, 59,60,61 i.e. 

= h {392 + 9 ' 2) " °- 044 ' (8 - 17) 



where g and g' are the SU(2)^ and U(l)y gauge couplings at the running scale 
Mz- The isosinglet heavy neutrinos also acquire thermal masses through Yukawa 
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Fig. 8. Lepton asymmetries for selected heavy Majorana neutrino scenarios. 



interactions 59 , i.e. 

m% (T) - m% (0) 1 
N ' T2 m = -^\hu\ 2 . (8.18) 

Such a T-dependcnt mass shift is small and comparable to the TV, widths at T w 
mjv 4 - Therefore, it is easy to see that the condition for resonant CP violation 
through mixing in Eq. (6.18) is qualitatively satisfied. Finally, the Higgs field also 
receives appreciable thermal contributions. The authors of Ref. 61 have computed 
the one-loop Higgs thermal mass, and found that M$(T)/T ~ 0.6 for values of the 
Higgs-boson mass favoured by LEP2, i.e. Mh < 200. In this range of Higgs masses, 
the thermal widths Tjy^T) will be reduced with respect to 1^(0) by a factor of 
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2 or 3 due to sizeable phase-space corrections. Nevertheless, the influence on the 
resonant phenomenon of CP violation through mixing is not dramatic when the 
latter effects arc included, and therefore large leptonic CP asymmetries are still 
conceivable. 
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9. Low-energy phenomenology of heavy Majorana neutrinos 

Whether heavy Majorana neutrinos can lead to interesting phenomenology at col- 
lider and lower energies is an issue that strongly depends on the out-of-equilibrium 
constraint given by Eq. (8.3). If this constraint is applied to all lepton families, heavy 
Majorana neutrinos have very little impact on collider phenomenology. However, 
this very last statement is rather model dependent. One can imagine, for example, 
a scenario in which AL e -violating operators exist and are out-of-equilibrium, and 
the fi and r sectors do not communicate any interaction to the electron sector, i.e. 
A(L e — L^) = A(L e — L T ) = 0. Since sphalerons conserve the individual quan- 
tum number B/3 — L e (see also Section 2), the observed baryonic asymmetry can 
be preserved in an excess of L e , independently of whether AL M - and/or AL T -non- 
conserving operators are in thermal equilibrium or not. As we will discuss below, 
such scenarios with strong mixing in the muon and tau sectors only can give rise to 
a variety of new-physics phenomena in a strength that can be probed in laboratory 
experiments. 

9.1. Lepton- flavour and/or number processes 

Heavy Majorana neutrinos with masses in the range 0.2 - 1 TeV may be produced 
directly at high-energy ee, 62 ep, 63 and pp colliders, 64 whose subsequent decays can 
give rise to distinct like-sign dilcpton signals. If heavy Majorana neutrinos are 
not accessible at high-energy colliders, they can still induce lepton-flavour- violating 
decays of the Z boson, 65 ' 66 the Higgs particle, 67 and the r and p leptons. 68 ' 69 
As we will see, non-decoupling quantum effects due to potentially large SU(2)l- 
breaking masses play a key role in these flavour-changing-neutral-current (FCNC) 
phenomena. 68 Heavy Majorana neutrinos may cause breaking of universality in lep- 
tonic diagonal Z-boson 70 and ir decays 49 or influence 71 the size of the electroweak 
oblique parameters S, T and U. 72 In fact, there exist many observables summa- 
rized in Ref. 73 to which heavy Majorana neutrinos may have sizeable contributions. 
These observables include r-polarization asymmetries, neutrino-counting experi- 
ments at the CERN Large Electron Positron Collider (LEP1) or at the Stanford 
Linear Accelerator (SLC), etc. 

In the following we shall show that high SU(2) ^-breaking masses in a class of neu- 
trino models can lead to large FCNC effects. Because these effects are not correlated 
with light neutrino masses, one can overcome the see-saw suppression relations that 
usually accompany such new-physics phenomena. 74 Let us consider a two-generation 
model of the kind. The model is similar to the one discussed in Section 3. It has 
two isosinglet neutrinos v' R and S' L . In the weak basis {{v^l) , { v tl) C , {S' L ) C , v'^) 
the neutrino mass matrix then takes the form 









m 2 M 



mi \ 
M 



M v = 






\ mi 



(9.1) 
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Diagonalization of M v yields two zero eigenvalues, which would correspond to mass- 
less fx and t neutrinos. If fi ^ and fi/M <C 1, they receive small radiative masses 
at high orders. 75 The other two states are very heavy, of order M±^t. In contrast to 
the usual seesaw scenario, the ratios rrii/M and m^/M remain fully unconstrained. 
Global analyses of low-energy data 73 restrict their values to mi/M, m^jM ~ 0.1. 
For reader's convenience, we define the parameters 



(ST) 2 - § , (#) a - ^ , (9-2) 

The newly introduced parameters quantify neutrino mixings between the light and 
heavy Majorana states. They also parameterize the deviation of the modified Wlvi 
coupling from the SM one. An extensive discussion is given in Ref. 73 ' 68 Including 
rcnormalization effects into the definition of light-heavy neutrino mixings, 49 one 
may tolerate the following upper limits 

(s^) 2 < 0.010, (s£) 2 < 0.035, and (s£) 2 < 1 x 10~ 8 . (9.3) 

The last limit comes from the requirement that only the electron family is respon- 
sible for baryogenesis. Of course, electron and muon families may interchange their 
roles in Eq. (9.3). 

Heavy Majorana neutrinos can induce sizeable FCNC decays of the type Z — > 
/it, r — > fi~e~e + or r — > fi~ fi~ /i + through quantum corrections presented in 
Figs. 9 and 10. Thus, the matrix element relevant for the generic decay r(p T ) — > 
l{pi)h{pi)h{p2) acquires contributions from 7- and Z-mediated graphs as well as 
from graphs with box diagrams. The respective transition elements are given by 



a 2 s 2 



<±M W 



*T 7 (r - Ihh) = S^W^l'X n F; l (j, - ^JL)(1 - 75) 



q 



-iG^a^{m T (l + 75) + m;(l - 75))] u T , (9.4) 

2 

iT z (t -> lhl 2 ) = y^r 8 hl2 Fl l u n ^{l - 75 K u h >f(l - As 2 w - 75 K , (9.5) 

2 

zT box (r -> J/Ja) = u l7 „(l - lb )u T u h ^(l - 75K , (9.6) 

where q = Pl + p 2 , s 2 w = 1 - M^/M 2 , and i^, G^, F^, F^ 1 ' 2 are certain 
composite form factors whose analytic form is given in Ref. 68 Nevertheless, it is 
useful to examine the asymptotic behaviour of the composite form factors for large 
values of \ Nl = m 2 Ni /M' 2 v and Ajv 2 = rr? N JM' 1 v in the two-generation model under 
discussion. For simplicity, we consider Ajv x ~ Ajv 2 ~ Ajv ~> 1- In this limit, we find 



F; 1 - --^^lnAw, (9.7) 



G;' - JaM, (9.8) 
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Fig. 10. Fcynman graphs pertaining to the decay r — ► I'lifo- 
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3 

F$ - -\sts v l\n\ N - Ists^^fXN, (9.9) 

i=l 

F T B U ol h - -(s-LS-26 hl2 +s^s^S U2 ) + \s v £s"ls v ts v t X N ■ (9.10) 

If all light-heavy neutrino mixings s£' are held fixed to a constant, the one-loop 
functions Gl^ F T Z \ and F^J 2 in Eqs. (9.7)-(9.10) do not vanish in the heavy 
neutrino limit, and therefore seem to violate the decoupling theorem due to Ap- 
pclquist and Carazzone. 76 However, it is known that the decoupling theorem does 
not apply to theories based on the spontaneous or dynamical breaking of gauge 
symmetries. Since we hold s2 fixed but increase the heavy-neutrino masses, the vi- 
olation of the decoupling theorem originates from the SU(2) ^-breaking Dirac mass 
terms m\ and TO2. 75 ' 67 The expected decoupling of the isosinglet mass M 77 can also 
be seen in Eqs. (9.7)-(9.10). This time we keep the Dirac masses mi and m 2 fixed 
and increase M. Taking Eq. (9.2) into account, it is then easy to show that all 
composite form factors vanish for large values of M « ttin 1 , "mN 2 ■ Consequently, 
there is a competitive loop effect of two scales, namely Dirac versus Majorana 
scale. High Dirac masses lead to non-decoupling loop effects while large Majorana 
masses give rise to a screening and reduce the strength of the effective FCNC cou- 
pling. Nevertheless, extensive analyses have shown that a non-decoupling 'window' 
confined by the two mass scales exists, within which FCNC and other new-physics 
phenomena come out to be rather large at a level that may be probed in next-round 
experiments. 67 ' 68 ' 69 

As examples, we calculate neutrinoless tau-lepton decays and flavour-violating 
Z-boson decays. Taking the dominant non-decoupling parts of the composite form 
factors into account, we arrive at the simple expressions for the branching ratios 



4 4 
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B(T ^ 6 C ) * 24576^ IV 
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+ 2(1-24)( s -) 2 ^( s -) 2 +84[E( s l) 2 ] 2 }, (9-11) 
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+ 2(1-24)(^) 2 ^( s -) 2 + 124[E( s l) 2 ] }> ( 9 - 12 ) 
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Fig. 11. Typical two-loop diagram contributing to the EDM of electron. 



* ^^(^(^[B^f , (9.13) 

where T T = 2.16x 10~ 12 GeV and Tz = 2.49 GeV are respectively the total widths of 
the t lepton and the Z boson known from experiment, 78 and T^^S) — F^f 1 /2. The 
complete analytic results of the branching ratios in Eqs. (9.11)— (9.13) are presented 
in Rcf. 68 

To give an estimate of the size of the FCNC effects, we take the maximally 
allowed values (s^) 2 = 0.035 and (s^) 2 = 0.010 given by Eq. (9.3). These light- 
heavy neutrino mixings lead to the branching ratios 

B(t~ fj, + ) ~ 2xl0~ 6 , B(t~ -» fi-e-e + ) ~ lxl0~ 6 , 

B(Z^fir) ~ 1.1 x 10~ 6 . (9.14) 

In Eq. (9.14), the upper limits are estimated by using the heavy-neutrino mass 
tojv — 4 TeV, which results from the requirement that perturbative unitarity be 
valid. The theoretical predictions of the branching ratios must be contrasted with 
the present experimental upper limits on these decays 78 

B(t~ -» B(t~ -> ^ere + ) < 1.4 x 10" 5 , 

B{Z^ht) < 1.3 xlO" 5 , (9.15) 

at the 90% confidence level. Future high-luminosity colliders and higher-precision 
experiments are capable of improving the above upper limits by one order of mag- 
nitude and so probe possible new-physics effects due to heavy Majorana neutrinos. 
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9.2. Electric dipole moment of the electron 

In general, CP-violating new-physics interactions may give rise to a large contribu- 
tion to the EDM of the electron. This results in an interaction in the Lagrangian 
of the form 



The experimental upper bound on electron EDM is very strict: (d e /e) < 10~ 26 cm. 78 
This bound is very crucial, as heavy Majorana neutrinos can induce an EDM of the 
electron at two loops. 79 A typical diagram is shown in Fig. 11. A simple estimate 
of this contribution based on a naive dimensional analysis for m^, S> M\y 

gives 21 



In Eq. (9.17) the factor depending on mjv 4 is always much smaller than unity. 
Clearly the above EDM limit could be important for \hu\ 3> 0.1 and/or ultra- 
heavy Majorana neutrinos with > 10 11 TeV. In this prediction, one should 
bear in mind that stability of the Higgs potential under radiative corrections requires 
\hu\ = 0(1). 80 Nevertheless, the EDM contribution is several orders of magnitude 

below the experimental bound for < 10~ 3 and/or \hi 2 \ <~ 10~ 2 . In this 

context, it is interesting to notice that lcptogcncsis models with nearly degenerate 
heavy Majorana neutrinos can naturally evade possible EDM constraints. 

10. Conclusions 

We have reviewed many recent developments that have been taking place in the 
scenario of baryogenesis through leptogenesis, and discussed the implications that 
heavy Majorana neutrinos may have for laboratory experiments. In the standard 
leptogenesis scenario, 15 L-violating decays of heavy Majorana neutrinos, which are 
out of thermal equilibrium, produce an excess in L that is converted into the ob- 
served BAU through B + L-violating interactions mediated by sphalerons. We have 
paid more attention to different kinds of mechanisms of CP violation involved in 
the Ni decays. One has two generic types: (a) CP violation originates from the in- 
terference between a tree-level graph and the absorptive part of the one-loop vertex 
(e'-type CP violation) and (b) CP violation comes from the interference between a 
tree-level graph and the absorptive part of the one-loop Ni self-energy (e-type CP 
violation). 

Recently, there has been renewed interest in £-type CP violation in models with 
mixed heavy Majorana neutrinos. If the masses of two heavy neutrinos become 
degenerate, then finite-order perturbation theory does no longer apply, and various 
methods have been invoked to cope with this problem in the literature. 20 Here, 
we have discussed the whole issue based on an effective resummation approach to 
unstable particle mixing. 21 One then finds that e-type CP violation is resonantly 




(9.16) 




(9.17) 
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enhanced if the mass splitting of the heavy Majorana neutrinos is comparable to 
their widths (cf. Eq. (6.18)), and if the parameter 5cp defined in Eq. (6.19) has a 
value close to 1. These two conditions turn out to be necessary and sufficient for 
resonant CP violation of order unity. 21 As a consequence, the scale of leptogenesis 
may be lowered up to TeV energies. In fact, Eg-motivated scenarios with nearly de- 
generate heavy Majorana neutrinos of order 1 TeV and universal Yukawa couplings 
can still be responsible for the BAU. This last observation receives firm support 
after solving numerically the BE's. 21 In this kinematic range, the e'-type contribu- 
tions are extremely suppressed. Also, finite-temperature effects on masses of heavy 
neutrinos and on decay widths do not spoil the above conditions for resonant CP 
asymmetries. Finally, constraints due to electron EDM are still too weak to play a 
role in leptogenesis. 

The fact that the isosinglet-neutrino scale can be lowered to TeV energies has 
a number of virtues. If one has to appeal to (local) supersymmetry in order to 
maintain the flatness of the inflaton potential, one then has to worry about the 
cosmological consequences of the gravitino during the nucleosynthesis epoch. Since 
the weakly interacting gravitinos are at most of the order of a few TeV, their slow 
decay rate will lead to an overproduction of D and 3 He unless the number density- 
to-entropy ratio at the time of reheating after inflation is less than about 1CP 10 . 5 
This leads to quite low reheating temperatures Trh ~ 10 10 GeV, after which the 
radiation-dominated era starts and baryogenesis or leptogenesis can in principle 
occur. The latter causes a major problem for GUT-scale baryogenesis and GUT's 
as well, especially if Trh ~ Mqtjt- In this context, it is important to remark that 
supersymmetric extensions of models with isosinglet neutrinos in the multi-TcV 
range as the ones discussed here can comfortably evade the known gravitino problem 
mentioned above. In such scenarios, the heavy inflatons with a mass of order 10 12 
GeV can now decay abundantly to the relatively lighter heavy Majorana neutrinos 
yielding equilibrated (incoherent) thermal distributions for the latter particles. 

Heavy Majorana neutrinos may also induce sizeable FCNC effects such as Z ^ 
Tfi and t — > /!/!/! at a level that can be probed in near-future experiments. Non- 
decoupling loop effects of high SU(2)i-breaking masses play a significant role in 
increasing drastically the strength of these new-physics phenomena. 75 ' 67 ' 68 ' 69 How- 
ever, these heavy Majorana neutrinos cannot account for the BAU at the same 
time. Depending on the model, they can coexist with the heavy Majorana neutri- 
nos responsible for leptogenesis without destroying baryogenesis. At the LHC, 64 
the viability of such models can be tested directly by looking for like-sign dilepton 
signals. Such signals will then strongly point towards the scenario of baryogenesis 
through leptogenesis as the underlying mechanism for understanding the baryonic 
asymmetry in nature. 
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Appendix A 

We shall now list useful analytic expressions for one-loop self-energies of the Higgs 
and fermion fields as well as for one-loop vertex couplings x + ZA^, x°viNi and HviNi. 
We present relations between wave-function CT's and unrenormalized self-energies 
in the OS renormalization scheme. The analytic results are expressed in terms of 
standard loop integrals presented in Ref. 81 Instead, we adopt the signature for the 
Minkowskian metric = diag(l, — 1, — 1, —1). 

The Feynman rules pertaining to the minimal model may be read off from the 
Lagrangian (3.7). We first give the the Higgs self-energies x~X~, X°X° an d HH, 
shown in Fig. 2(d)-(f) 

n x - x -(p 2 ) = n x o x o( P 2 ) = n HH (p 2 ) 

= EE TT [<#o(p 2 ,r<,0) + p 2 B 1 (p 2 ,m%,0)].(A.l) 

1=1 i=l 7F 

From Eq. (A.l), the universality of the divergent parts of the wave functions 8Z X - , 
SZ x o and SZh is evident, since 5Z$ V = —RcH'£ lv (0) for all field components of the 
Higgs doublet <f>. 

Figures 2(g), (h) and (j) show the individual contributions to the one-loop 
fcrmionic transitions, V — ► I, vy — » v\ and Nj — > Ni, respectively. Explicit cal- 
culation of these self-energy transitions gives 

i—1 

+ B 1 {p 2 ,m 2 Ni ,M 2 H )^ + m Ni (h* u h* Vi P R + h H h Vi P L ) 

x (^(p 2 , 774,0) - Bidi 2 ,^,^))], (A.3) 

^n %N] (f) = - E 1^2 ( h ti h U * p R + huhlj *Pl) ( (P 2 ' 0, 0) 
1=1 11 

+ l -B 1 { P \Q,Ml)). (A.4) 

Note that the light-neutrino self-energies in Eq. (A.3) contain non-zero masses in 
the limit /5 — s- if My 7^ 0. Therefore, at T = 0, small radiative neutrino masses 
can be generated. 75 However, these contributions are suppressed by small Yukawa 
coupling, and therefore do not invalidate any experimental or cosmological limit. 

Before we express the wave-function renormalization factors in terms of un- 
renormalized self-energies, we first notice that the one-loop fermionic transitions 
fj — > fi (with fi = I, i>i, Ni) given by Eqs. (A.2)— (A.4) have the generic form 

E y (jQ = £f> 2 ) $P L + S«( P 2 ) ffP R + Z%(j?)P L + Z%*(p 2 )P R , (A.5) 



Heavy Majorana neutrinos and baryogenesis 47 



where only dispersive parts are considered. If the transitions involve Majorana 
fermions only, one then has the additional properties T<fj(p 2 ) — Y>f*{p 2 ) and T,fj (p 2 ) - 
E^f (p 2 ). Following Rcf. 49 , the wave-function CT's are given by 



SZL 



and, for i ^ j, 



+ 



-E^(m 2 ) 2m 2 Ef i '(m 2 ) m^Sf (^) +Sf (m?) 
1 

2rrii 



(A.6) 



5Zl 



m\ - m 2 



'"JE';!'";) + m^E^m*) + m^{m^ + m^E^O 



Mi 



(A.7) 

The wave-function renormalization of charged lcptons may be obtained by Eqs. 
(A.6) and (A.7), if all terms depending on E^ (p 2 ) and its derivative are neglected. 
At this point, it is important to remark that there will be additional contributions 
to our wave- function CT's of the Higgs and fermion fields coming from loops related 
to SM gauge particles. As we have seen in Section 8, the same kind of loops can 
induce non-zero thermal masses to the Higgs and fermion particles. 59,60 These new 
contributions are universal and therefore pose no problem to the Yukawa-coupling 
renormalization discussed in Section 4. 

Finally, one-loop corrections to the vertices x^^A^, x° v i^i an d Hv\Ni shown 
in Figs. 2(a)-(c) have been calculated in Ref. 21 by keeping the complete functional 
dependence on the Higgs-boson mass Mh- Their analytic expressions are given by 



iV 



-iUlP R U Ni 2^}^ 167r 2 h l'i h l'3 h lj 
1<=1 j=l 

x C (0,0,m 2 ,0,mi.,0) + C 12 (0, 0, m 2 , 0, m 2 . , 0) 



(A.8) 



x C (0,0,m 2 Ni ,0,m%.,0) + C Q (0,0,m 2 Ni , M 2 ,m%.,0) 
+C 12 (0,0,m 2 Ni ,0,m%.,0) + C 12 (0,0,m 2 Ni , M 2 ,m%.,0) 

x m 2 Ni (c 12 (0, 0,m 2 Ni ,0,m 2 N .,0) - Ci 2 (0,0,m^,M^,m^.,0)) 
-M 2 C (0,0,m 2 Ni ,M 2 ,m 2 N .,0)\ }, (A.9) 



-iv Hvi N, = iuiP R u Ni I 



m Nt m Nj 
32V2tt 2 



h* Vi h V jhij 



V = lj=l 

C (0,0,m 2 ,0,m 2 .,0) + C (0, 0, m% , M^, m% , 0) 
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+C 12 {0,0,m 2 Ni ,0,m 2 N .,0) + C 12 (0,0,m^,M|,m^.,0) 
H ?= — hp ihp aHu 



x 
- M 



m 2 Ni (c 12 (0,0,m 2 Ni> 0,m 2 N .,0) - Ci 2 (0,0,m^,M^,m^,0)) 
\I 2 C (0,0,m 2 Ni ,M 2 ,m 2 N .,0)\ }. (A.10) 

From Eqs. (A. 9) and (A. 10) it is easy to see that the L-conserving part of the 
couplings x°^;iVi and HviNi is UV finite and vanishes identically for Mh — > 0. 



